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scale is essentially given by the atom rest frame redistribution asymptotic. The pure Doppler 
and combined natural, Doppler and collisional broadening are characterized by both small 
and large scale superdiffusive Levy flight behaviors while the combined natural and Doppler 
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in the wings of the spectral distribution is incompatible with the breakdown of superdiffusion 
in combined natural and Doppler broadening conditions. 
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INTRODUCTION 



The study of transport phenomena [1] and in particular of diffusion processes {2] has 
developed considerably in the course of the last decades. One of the foundations of this 
area is the well-known diffusion equation, consequence of the central limit theorem [3j, that 
describes the random motion of a particle in an isotropic homogeneous three-dimensional 
space. The Probability Density Function (PDF) of the diffusing species being at a certain 

y high number of jumps and for an initial location 
5] 



position r at time t, after a sufficient 



at r = 0, is given by the propagator 



1 / 3r 2 12 



with the root mean square displacement, 



(r 2 (t)) = 6Dt, (2) 
where D is the diffusion coefficient of the diffusing species. 

By the end of the last century it was recognized that the simple behavior in Eqs. (H|) 
and ([2]) was not completely general. The two most widely known deviations from the standard 
diffusive behavior are Levy flights [g, 3], coined and popularized by Benoit Maldelbrot, and 



Richardson's law of turbulent diffusion [8]. Recently, anomalous diffusion equations were 
also found to apply to in the econometric modeling and/or prediction [9], in laser cooling 
of atoms [llJ], in the foraging behavior of some animals [ill ], and to the the scaling laws 
ruling human travel activities (critically connected to the geographical spreading of infectious 



diseases) [12]. 

Situations deviating from the classical diffusion case can heuristically be described by a 
modified power law giving the scaling of the mean squared displacement as 



<r 2 (t))oc^ (3) 

Anomalous diffusion results in the departure of the 7 factor from unity in the last equation, 
owing to the non applicability of the classical central limit theorem, as a result of the presence 
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of broad distributions or long-range correlations and non-local effects [j], [3] . Broad spatial 
jump or waiting time distributions lead to a non-Gaussian and possibly a non-Markovian 



time and spatial evolution of the system 
be collectively described as strange kinetics 
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giving rise to anomalous behaviors that can 



Equation ([3]) encompasses both the linear dependence of the mean squared displacement 
with time (7 = 1) characteristic of the standard Brownian motion as well as nonlinear de- 
pendencies, slower (7 < 1) and faster (7 > 1) than the classical case. The below Browian 
motion regime is usually named sub-diffusion or dispersive motion while the above range is 
known as hyper or superdiffusion. In order to generalize the classical treatment to consider 
anomalous diffusion, the Continuous Time Random Walk (CTRW) model is often considered. 
The random trajectory is viewed as the result of pairwise stochastically independent spatial 
and temporal increment events. When both the variance of the spatial steps and the expecta- 
tion value of the temporal increments are finite, CTRW is equivalent to Brownian motion on 
large spatio-temporal scales [5j and yields ordinary diffusion. Nevertheless, for the anomalous 
diffusion, a bifractional diffusion equation can be derived for the dynamics of W (r,t). This 
can be solved using the methods of fractional calculus to give 

where is a one-sided Levy stable law of index fi [13]. This equation implies that the 
mean squared displacement scales as 

p(r) ~ r - {l+ ^. (5) 

The particle's motion can be classified either as a Levy flight, in the case that the time- 
of-flight of each jump in the trajectory is negligible, or as a Levy walk, if the finite time to 
complete each jump must be taken into account P. w\. 

In this work we are interested in the stochastic theory of atomic resonance radiation mi- 
gration under which the trajectories are of the supperdiffusive type. The radiation transport 
is usually described in terms of a master equation, which has an integro-differential form 
taking into account non- 



equation, named after seminal contributions in the 1940s [r]. fis| . It was shown, first using 
ad hoc arguments [19|] and later demonstrated, that the motion of photons is super diffusive. 
This was initially demonstrated for Doppler and Lorentz spectral line shapes and recently 
generalized for any line shape [151 ]. The work was done in the limit of Complete Frequency 
Redistribution (CFR), a situation in which the collision rate is sufficiently high when com- 
pared to the lifetime of the excited state to destroy any correlation between photon absorption 
and (re)emission events. This contribution aims to extend the previous work and discuss the 
influence of Partial Frequency Redistribution (PFR) effects in the stochastic description of 
the resonance radiation trajectories. We will show that the overall asymptotic behavior of the 
trajectories is dictated by the redistribution in the rest frame of the atom, even if in the lab 
rest frame memory effects persist in a finer scale. Specifically, we will show that the Doppler 
or Lorentz CFR superdiffusive asymptotics are recovered for PFR if there is only Doppler or 
Doppler plus natural redistribution in the atom frame. On the contrary, we show for the first 
time, that the complete coherence in the atom's rest frame will manifest itself as a breakdown 
of the superdiffusive behavior in all cases and not only for the wing frequencies events. 

In Sec. II the Partial Frequency Redistribution (PFR) effects are introduced along with 
the notation to be used. In Sec. Ill the single line Opacity Distribution Function (ODF) 
under PFR is defined. This is subsequently used to derive the the large scale asymptotics of the 
jump length distribution in Sec. IV and of the opacity distribution in Sec. V. Section VI makes 
a connection between the superdiffusive character and the low opacity asymptotic behavior 
of the ODF. In Sec. VII numerical results are presented for the ODFs and the conditional 
jump length probability distribution functions, and the influence of frequency redistribution 
in both a small scale (controlled by PFR effects) as well as in a large scale (corresponding 
to either a CFR superdiffusive or, in alternative, a diffusive) asymptotic are discussed. The 
main conclusions are summarized at the end (Sec. VIII). 

II. PARTIAL FREQUENCY REDISTRIBUTION 

Resonance radiation trapping can be envisaged as a random flight with a jump size distri- 
bution dependent upon the spectral line shape. If all the scattering (absorption followed by 
reemission) events are coherent in frequency, the radiation transport is described by a diffu- 
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sion equation with a diffusion coefficient dependent on the mean free path of the radiation. 
However, reabsorption-reemission events are inelastic, and a photon frequency redistribution 
in the lab reference frame exists preventing the notion of a mean free path. Therefore, the 
topology of the excitation random trajectory has to be related with the frequency redistribu- 
tion of the emitted photons, since the jump length Probability Distribution Function (PDF) 
must consider all possible optical emission frequencies. This can be done as 
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»+oo 

p(r) = I © (x) p (r\x) dx 

-oo 
-oo 

6 (x) $ (x) e-* (x)r dx, (6) 

-oo 

where x is a frequency difference to the center-of-line frequency x = , properly nor- 
malized by dividing by the corresponding width parameter (given for the Doppler and Lorentz 
spectral profiles as Woop = v mp ^ and Wlot = r/47r, respectively, with v mp the most probable 
thermal velocity and T the radiative rate constant), G (x) is the emission spectrum lineshape, 
p (r\x) the photon exponential (Beer-Lambert) jump distribution, and (x) is the normalized 
absorption lineshape of the resonance line( so that 3> (x) dx). The dimensionless distance 
r = defines an opacity or optical density scale from the center-of-line optical depth ko 
and the center-of-line absorption coefficient (0). For a sufficiently high number of collisions 
during the lifetime of excited atoms, the emitted photon has a frequency completely uncor- 
rected with the absorption frequency and so the radiation is completely redistributed over 
the entire spectral line. This corresponds to the Complete Frequency Redistribution (CFR) 
case, for which the emission and absorption line shapes coincide and the jump length PDF is 
independent of the past history. For strong resonance lines with short lifetimes the frequency 
redistribution is partial since only a few collisions occur before emission, and consequently 
the stochastic nature of the distribution must be considered. Partial Frequency Redistribu- 
tion (PFR) effects are important in an astrophysical context, as is the case of the scattering 



of Lya radiation in optically thick nebulae 



20J and the Ca I resonance line in the solar 



spectrum [211 ] . as well as in lab scale atomic vapors, notably the 185 nm mercury line in 
low-pressure lighting discharges (which can account for as many as 10% of the overall flux 
in a typical T12 fluorescence lamp) [3] and the 147 and 129 nm Xe VUV radiation used 
in Plasma Display Pannels (PDPs) applications [2^ . For both mercury and xenon, the 



5 



natural lifetimes of the corresponding excited states are less than about 3 ns and therefore 
only at very high densities is the collision rate high enough to approach CFR conditions. 
For a two-level model in the CFR limit, both absorption and emission spectra lineshapes 



or 



can be described by Doppler - $_o op (x) = ^e x -, Lorentz - &L or (x) = -j^p 
Voigt - &voigt (x) = ^§j2 I-oo a 2+( x u ^ du ~ spectral distributions, in which x is a normal- 
ized (with the Doppler width) difference to the center of line frequency and a is the Voigt 
characteristic width, a ratio of the Lorentz to the Doppler widths defined previously. The 
study of PFR was pioneered in the 1940s in an astrophysical context and we will use the 
notation introduced by Hummer in the early 1960s [25|. PFR is described by a redistribution 
R function which is the joint probability of absorption of a x' photon and reemission of a 
x frequency photon. This redistribution function is derived for a thermal vapor, assuming a 
Maxwell-Boltzmann velocity distribution for ground state atoms and further imposing that 
the atoms velocity is unchanged upon excitation (assumption not valid for ultracold vapors). 

We will further consider unpolarized radiation (scalar redistribution function), use the an- 
gle averaged redistribution for an isotropic angular distribution reemission and discuss the 
three most important cases for a ground state absorbing state: (i) pure Doppler broaden- 



ing (Ri), (h) combined natural and Dopp 
natural and collisional broadening (Rui) 



er broadening (Rn) and (iii) combined Doppler, 



16, 



20, 



26, 



271 ] . For the Rj case, both levels are 



infinitely sharp and there is pure coherent scattering in the Atom's Rest Frame (ARF). This 
is unrealistic because both the lower and upper levels of the optical transition are broadened 
to some extent by radiation damping and/or collisions. It is nevertheless useful since it allows 
the separation of the Doppler shift from the other broadening mechanisms. The physical 
picture for Rn is a line with an infinitely sharp lower level and an upper level broadened only 
by radiation damping. In the ARF the absorption is Lorentzian and scattering is completely 
(frequency) coherent. This type of scattering is important for resonance lines in low density 
media where the collisions do not perturb the excited levels. Rni applies whenever the up- 
per state is broadened by radiation and collision damping in the limit where collisions are 
frequent enough to cause complete frequency redistribution in the ARF. Both absorption and 
reemission profiles are Lorentzian in the atom's frame with a total width equal to the sum 
of the radiative plus collisional (uncorrelated) characteristic widths. It is well known that 
complete redistribution in the ARF does not imply the same quantitative behavior in the lab 
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frame used to describe radiation transport. Under the 
lab frame joint redistribution functions are given by 



previously mentioned assumptions, the 
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28] 



Ri (x'; x) = - erfc (x) 



(7) 



R H (x';x) 



IT' 



3/2 



+oo 



l\x-x'\ 



x < arctan 



x + u 



a 



arctan 



x — u 



du, 



(8) 



and 



Rin (x; x) 



1 



r 5/2 



+oo 



x { arctan 



x + u 



x < arctan 



a 

x + u 
a 



arctan 



x — u 



(9) 





x — u 


— arctan 






a 



du, 



where a is the total width, erfc is the complementary error function, and x and x stand 
for x = min , \x\) and x = max (|x'| , |x|). 

In Eq. ([6]) we need to use the (re)emission PDF. For radiation transport in thermal va- 
pors we can assume a CFR absorption profile which corresponds to the ensemble average 
over the Maxwellian velocity distribution. The reemission distribution comes from the joint 
redistribution probability and, according to Bayes rule, is given by 



9 (x\x') 



R (x'; x) 
$(x') 



(10) 



and it is this conditional (re)emission PDF that should be used. <3? (x 1 ) is the absorp- 
tion spectrum PDF which should be either the pure Doppler distribution or the Voigt 
profile (Rn and Rin). 

Figures Q] to [3] show the spectral distributions calculated from the previous equations. Both 
Ri and Rni distributions are symmetric, while Rn is asymmetric. In all complete 
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FIG. 1: Conditional reemission Probability Distribution Functions (PDF) of Partial Redistribution 
function Rj for several previous absorption frequencies x'. The corresponding Doppler Complete 
Redistribution is also shown. 

redistribution in the ARF does not imply CFR in the lab frame, the degree of redistribution 
in the lab frame being always bigger for core radiation and for the Ri and Rni cases. The 
complete coherence in the atom's rest frame for the Rn model is not complete in the lab frame. 
However, the coherence degree is more retained for wing photons and higher a values (bigger 
importance of natural coherent broadening compared to Doppler). On the other hand, for the 
Rin distribution, the increase of the a value implies higher redistribution in the core (bigger 
importance of collisional redistributed broadening). We will discuss now the influence of the 
redistribution between absorption and reemission in the asymptotic super diffusive character 
of photon migration under PFR conditions. For this purpose we extend the concept of a single 
line opacity distribution function for incomplete redistribution. 




III. LINE OPACITY DISTRIBUTION FUNCTION UNDER PARTIAL 

REDISTRIBUTION 

To study the super diffusive behavior of resonance radiation trapping it is convenient to dis- 
cuss the influence of the line shape in the Opacity (probability) Distribution Function (ODF). 
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FIG. 2: Conditional reemission Probability Distribution Functions (PDF) of Partial Redistribution 
function Rm for several previous absorption frequencies x' and Voigt width parameters a. The 
corresponding Voigt Complete Redistribution is also shown. 

For that we will use a line opacity scale, defined from the normalized absorption line shape 
as k = (x). This line opacity scale will thus have no reference to the actual system size 
or number density and will be useful only for infinite media. This is the one chosen to 
discuss the ODF. The monochromatic line opacity along a given pathlength I, for homoge- 
neously distributed particles, is k (x) = uctqI^ (x) /$ (0) = kqQ (x) /$ (0) = <J> (x) r and is 
therefore proportional to the number density n (or the pressure) and to the center-of-line 
opacity ko = na^l. ctq is the center-of-line absorption cross section as usual. The discus- 
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FIG. 3: Conditional reemission Probability Distribution Functions (PDF) of Partial Redistribution 
function Rn for several previous absorption frequencies x' and Voigt width parameters a. The corre- 
sponding Voigt Complete Redistribution is also shown. 

sion of the spectral line profile influence on the jump length distribution will use an overall 
(for the whole line) dimensionless opacity or optical density scale defined as r = so 
that r = k (x) dx. This is done since any given jump distance must reflect the whole of 
the line shape and not just a given monochromatic opacity. 

We will now generalize the basic procedure used previously 15[] for the complete frequency 



distribution to cover the partial frequency distribution (PFR) case. We will start with the 
jump length PDF by rewriting Eq. (|6|) making explicit the conditional dependence in the 
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absorption frequency: 



/+oo 
@(x\x')®(x)e-*( x >dx. (11) 
-oo 

Since photon migration by reabsorption is associated with exponential jump size PDFs, 
it is desirable to consider the general p (r\x') as a linear combination of exponential densities 
given by 



-oo 

-ft r 







p(r\x') = H(k\x') (ke- hr ) dk, (12) 



where H (k\x ! ) is the PDF of effective line opacities. The single line opacity at a given 
frequency, x, is given by k (x) = <3? (x). As was shown before for the CFR case, kH (k\x') is 
the inverse Laplace transform of p(r\x'). The inversion is analytical and can be performed 



using the real inversion form of the Laplace transform (with the parameter c = 0) 29f] , 



Y r+co 

k H (k\x) = — / {Re \p {iui)\ cos {koj) — Im \p (iu>)] sin (ku>)} duo, (13) 
71 Jo 

which gives 



i r+oo 

H (k\x) = - I 9 (x\x') $(x)6[k-& (x)] die. (14) 



oo 



This can be further simplified by making the change of variable y = $ (x) and decomposing 
the integration into positive and negative x frequencies, 

1 /■*((>) 

H (k\x') = ~iJ o [0+ (k) + Q- (*)] y *' (y) S[k- y ] dy, (15) 

since it is assumed that $ (x) is always nonnegative, and that the transformation has 
opposite signs for positive and negative frequencies. (y) is the inverse function of ^(x), 
\J/ (y) = (y). The integration in y only runs to the maximum value of the absorption 
spectrum (the center- of -line absorption coefficient, ^(0)) and 6+(k) and 9- (k) represent 
the value of the (conditional) probabilities for emission for positive and negative frequencies, 
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implicitly given from the opacity k through the absorption spectrum. For any particular 
opacity k, there are two (symmetrical) frequencies, +x and —x, that give rise to k defined 
by (x) = k (symmetric absorption). These frequencies do not necessarily have the same 
reemission probability. One therefore can define 9+ (k) = Q (x\x') and 9- (k) = O {—x\x') and 
finally rewrite Eq. (|15p as 

f -[M*0+M*)]*'(*0 if k<$(0) 
H [k\x J = < • (16) 

[ if jfe > * (0) 

This is the general workhorse that will be used in the following sections to discuss the 
influence of PFR in resonance radiation reabsorption and compare it with the CFR case fl^ . 
We will show how the redistribution of the emission into the wings of the absorption spectrum 
implies a superdiffusive behavior characterized by a very broad (infinite mean) jump size 
distribution. In addition, we are able to demonstrate that, if redistribution into the tails of 
the absorption is prevented due to PFR memory effects, the superdiffusive character breaks 
down, and this can be seen in a cut-off in the PDF of the opacity at small line opacity values 
that is translated into a corresponding cut-off of the higher jump lengths. 

Two of the PFR cases defined before (Rj and Rni) have symmetric reemission PDFs (see 
Figs. 1 and 2), which allow us to write, 



H (k\x') 



-29+ (k) (k) if k < $ (0) 
if Jfe > $ (0) 

For Complete Frequency Redistribution this gives 



(17) 



. -2Jfe¥'(Jfe) if k <$(0) 
H{k\x') = { , (18) 

if k > $ (0) 



which is simply the result obtained before [151 ] . 



IV. ASYMPTOTIC BEHAVIOR OF THE JUMP DISTRIBUTION 



Using Eq. (|16p . Eq. (|12p can be rewritten as 
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*(0) 



p (r\x') = - [9 + (k) + (jfe)] tf' (jfe) (jfe e~ fer ) dfc. 



(19) 



The asymptotic behavior is easily found by performing the change of variable y = kr, 



p (r\x') 



*(0)r 



which gives 



+ 



[ye 



') dy, 



(20) 



p (r\x') ~ r 2 



r J \r 



y'( U ~) {ye~y) dy. 



(21) 



For PFR conditions, we have to consider both the absorption as well as the conditional 
reemission asymptotics. For the absorption spectrum, if the line shape function can be given 



by a power-law asymptotic in the wings 



15, 



191 ]. then one has 



k(x) = 3>(z) ~ x Pa with (p a > 1), 



where a stands for absorption. 
The inverse function is 



^ (k) ~ k Pa 

k— >0 



which can be written, as 



T , , y\ r _j_ _l 

V|/ [ — ) r^i y Pa J- Pa , 

,r/ r^oo p a y 



(22) 



(23) 



(24) 



If the asymptotic of the conditional reemission is similar to that of absorption (but possibly 
with a different parameter value), then 



(x\x') ~ x Pe with (p e > 1) 



(25) 
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where e stands for emission. Therefore, both 9 + and (9_ scale as 



?+(*;) ~ 6L (k) ~ fc^ 



(26) 



Finally, Eq. (|21|) can be rewritten as 



P { r \ x ') 



-2 



+oo 



y\Pe/Pa 



' X -l/p ar X/p a 
PaV ' ' 



(ye y ) dy 



1(1 + — - —| r p « 

Pa Pa 



(27) 



and therefore, 



/ I A — - fl+^l 

p I T VX J ~ r Pa ^ Pa / 1 



(28) 



Comparing the previous equation with the Levy flight asymptotic in Eq. ([5]), gives 



1 , 

M = — (Pe — 1) 
Pa 



(29) 



If the asymptotic regime is the same for both absorption and PFR reemission, p = p a = p e 
and the jump size distribution has the asymptotic 



P { r \ x ') 



-(2-1) 



(30) 



which corresponds to the CFR form obtained before 15|]. So, the CFR asymptotic is 
recovered even in PFR conditions, as long as the asymptotic forms of both the absorption 
spectrum and the conditional reemission are the same. The Levy parameter is also the one 
obtained in CFR conditions, 



P 



(31) 
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V. ASYMPTOTIC BEHAVIOR OF THE OPACITY DISTRIBUTION 



We can obtain the same conclusions from the spectral and from the ODF PDFs asymp- 
totics. Indeed, the small opacity scaling law can be easily obtained from the previous equa- 
tions. From Eqs. (|22p and (j23[) . one obtains 



*' (k) 



— k p<^ 



k^O p a 

Using the reemission asymptotic of Eq. (|26p . the ODF small opacity limit is 



(32) 



H(k\x') — k ~ 



)]. 



(33) 



which, for absorption and reemission spectral distributions with the same asymptotic char- 
acter p = p a = p e , gives 



H(k\x') ~ —k p. 



(34) 



' fc^o p 

The last two equations show that the ODF asymptotic behavior can be easily identified 
from a log-log plot as 



log if (jfelx') 



k^O 



Pa 



+ 1 



Pe 
Pa 



and 



log A;, 



(35) 



log if (k\x') ~ — log A;. 
fc^O p 



(36) 



The cases of CFR Doppler and Lorentz resonance line radiation transfer give analytical 
ODF expressions and were considered in previous works [15I. lit]]. For the Doppler case, one 
has fj, ~ 1 and thus Eq. (|22p is only approximately valid and the dependence of the ODF 
over the opacity values is extremely weak. However, for Lorentz, the asymptotic in Eq. 
is exact with p = 2. The ODF is given by [If 
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^w-TCT- 1 )"""- (37> 

which has an exact asymptotic of the form given by Eq. (|36p . 

VI. SUPERDIFFUSIVE BEHAVIOR FOR RESONANCE LINE TRAPPING 

Although the super diffusive behavior is well established from the asymptotics derived in 
the previous two sections, a simpler more informative, alternative approach is possible. Equa- 
tion (fT2l) can be written in Laplace space as 

, . A kH(k\x') „ , . 

p(ux') = / y —^dk, (38) 

K ' Jo u + k y ' 

and since the moments of the jump size distribution can be recast from the mth derivatives 

of this Laplace transform as p( m ^ (u\x')\ u _ Q = (— l) m {r m ), one obtains 

/ r m\ =m] n W- 1 ) dk = m] / k -m\ (39) 

Jo K 

The integration only goes up to the upper limit of the center-of-line opacity. The first 
moment gives (r) = J* *^ H ^ x ) d/c which shows that, unless H (k\x') vanishes as k — ► 0, an 
infinite mean jump size is obtained. 

If one considers now the asymptotic of Eq. (I34p . one can conclude that (r) oc lim ttt- 
which is finite only if 1/p < 0. But, since p > 1 by construction, one can finally conclude that 
the mean jump size is infinite irrespective of the actual value of the p parameter. 



VII. NUMERICAL RESULTS 

Figures H] to [6] give the conditional ODFs for the three PFR redistribution mechanisms 
considered in this work alongside with the corresponding CFR Doppler and Lorentz limiting 
cases. The results were computed from Eqs. ()16p and ()17p . The CFR ODFs allow one 
to quantify very easily the cause of the super diffusive behavior for resonance line radiation 
trapping; for very small opacities, the ODF either decreases very slowly (Doppler) or even 
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increases in importance (Lorentz) with a decrease in opacity values. It does not vanish in 
the limit of k — > and therefore superdiffusion sets in (see previous section). From Eq. (|12p . 
one can see that the slower the decrease in the ODF for smaller opacity values, the more 
important the higher jump sizes. For Doppler, Lorentz (and Voigt) CFR trapping, the ODF 
variation for small k values is either increasing or so slowly decreasing as to render even the 
first moment of the jump size distribution infinite [15I. [ig|. The CFR Voigt case is a self-affine 
multifractal with two different characteristic Levy flight fj,'s parameters, a fj, ~ 1 corresponding 
to the Doppler-like core radiation and a fj, = 1/2 for Lorentz-like wing photons. The overall 
excitation random walk topology is controlled by the bigger Lorentz scale (if not truncated 
by the vapor confining boundaries) [l9 |. Therefore, in the present case of PFR effects in an 
infinite medium, and even for the Voigt Maxwellian absorption profile, we need only concern 
ourselves with the PFR deviation from the most extreme Lorentz CFR asymptotic. 

Figure H] shows that the superdiffusive character of Ri partial redistribuion is similar to the 
one for the corresponding Doppler CFR case. There is a transition between two regimes: (i) the 
central core constant redistribution (compare Fig. [1]), characterized by a faster than Lorentz 
increase of opacity PDF for small opacity values, and (ii) the wings Doppler-like redistribution, 
with a Doppler-like (albeit with a slower convergence) asymptotic. The transition between the 
two regimes depends upon the actual value of the previous absorption frequency x' (for x' = 4 
and 5, the transition still occurs but at opacity values smaller than the ones represented in 
Fig. i) . 

The Ri case is useful because it allows one to single out the pure Doppler effect on the 
redistribution functions. The more realistic cases are described by the Rn or Rni distri- 
butions. Figure [5] show the conditional ODFs for the Rni redistribution for representative 
values of the a parameter found in experimental conditions. The PFR random walk has a su- 
perdiffusive character with the same asymptotic of the CFR Lorentz. For the Rni combined 
Doppler, natural and collisional broadening, there is complete (Lorentzian) redistribution in 
the atom's rest frame, but not in the lab frame. One can distinguish mainly two regimes, one 
for very small line opacities and the other for opacities approaching the center-of-line value, 
the exact transition between the two depending upon the a parameter. For very small line 
opacities, the actual superdiffusive regime corresponds to the CFR Lorentz case. For higher 
line opacities, the redistribution is similar to Ri; an initial step increase in ODF PDF with 
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FIG. 4: Conditional Opacity Distribution Function (ODF) for Partial Redistribution function Rj for 
several previous absorption frequencies x' . The Complete Redistribution limiting Doppler and Lorentz 
cases are also shown. 

decreasing opacity, followed by a Doppler-like asymptotic (which eventually merges into the 
Lorentz higher scale asymptotic [l5]). The a parameter quantifies the relative importance 
of Lorentz (natural plus collisional in this case) over Doppler; the higher the a value the 
more important the Lorentz complete redistribution in the atom's rest frame and therefore 
the smaller the deviations from CFR and the onset of the CFR Lorentz asymptotic occurs 
earlier (ie., for higher opacity values). 

The conditional ODFs for the Rjj redistribution are shown in Fig. [6l The most remarkable 
difference from the previous cases resides in the low line opacity asymptotic. There is a very 
abrupt cut-off of the ODF PDF values for small opacities which scales approximately as a. 
This case corresponds to excited levels not perturbed by collisions during their lifetime and 
coherent reemission in the Atom's Rest Frame (ARF). This coherent reemission in the ARF 
does not give automatically coherence in the Lab Rest Frame (LRF) due to the Maxwell 
velocity distribution. However, after several Doppler widths (depending on a), the ARF 
coherency appears in the cut-off of the ODF. As for the previous Rni case, two regimes 
can be distinguished, with a transition between them that depends on a. For line opacities 
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FIG. 5: Conditional Opacity Distribution Function (ODF) for Partial Redistribution function Rju for 
several previous absorption frequencies x' and Voigt width parameters a. The Complete Redistribution 
limiting Doppler and Lorentz cases are also shown. 
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approaching the center-of-line value, $0) the redistribution is similar to the pure Doppler 
Ri case; an initial very steep (faster than Lorentzian and persisting to smaller line opacities 
the more off-center is the previous photon absorption) increase with decreasing opacity value 
followed by a CFR Doppler redistribution until the small opacities approaching the cut-off 
values (compare the CFR asymptotic with the ODF curves for the x' = 0, 1 and 2 cases). 
The cut-off of the ODF PDF for the smaller line opacities is very important since it shows 
the breakdown of the superdiffusion character for this Rn model. 

From the above results some general conclusions can be extracted: (i) the ARF frequency 
coherence does not immediately give rise to LRF coherence due to the ensemble Maxwellian 
velocity distribution averaging and (ii) the ARF redistribution (or non redistribution) char- 
acter will eventually show up as a corresponding LRF asymptotic for sufficiently small line 
opacities (that is to say, for optically thick enough vapors). The CFR in the ARF for the Rni 
will give rise to a Lorentz CFR (superdiffusive) asymptotic in infinite media even under par- 
tial redistribution while the complete coherence in the ARF for Rn will reveal itself in the 
LRF as a breakdown of the superdiffusive charater. The breakdown of the superdiffusion for 
thick vapors under Rn conditions is particularly important due to the widespread practice 
of reducing the actual complex behavior of the Rn redistribution into a linear superposition 
of complete frequency redistribution in the line core plus a complete coherence in the line 
wings. This practice is essentially motivated by the mathematical simplification that it allows 
and it can be traced back to the work of Jefferies and White (and it is therefore known 
as the J W- approximation). Although the several implementations of the JW-approximation 
differentiate the photons from the core from those from the wings, all ignore the possibility 
that even the most extreme wing photons can be redistributed in frequency, thus neglecting 
the so-called "difusion" in frequency. Moreover, the full redistribution function is skewed 
towards the line center and this is not taken into account. Figure [6] clearly shows another 
physical limitation of the J W- approximations which is generally not pointed out. For high 
opacity samples, the CFR assumption in core radiation is a severe approximation since CFR 
implies that radiation well into the wings can be emitted in a single scattering event, which 
can not occur (an upper limit exists) due to the breakdown of the Doppler asymptotic. This 



contributes to show that the JW-type approximations are only acceptable for sma 
parameters and for overall center-of-line opacities roughly smaller than 100 — 500 [r] . 
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FIG. 6: Conditional Opacity Distribution Function (ODF) for Partial Redistribution function Ru for 
several previous absorption frequencies x' and Voigt width parameters a. The Complete Redistribution 
limiting Doppler and Lorentz cases are also shown. 
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The asymptotic analysis of resonance radiation trapping started probably with Holstein 
seminal work [2] and received afterwards a relevant contribution from van Trigt's infinite 



opacity expansion analysis in the 1970s 



321 ] . Asymptotic analysis is a very useful tool to 



calculate the large-scale behavior, which, from the derived scaling laws, can provide physical 
insight into more complex and realistic conditions. This is particularly true for partial redis- 
tribution in the LRF. We clearly show that the large scale behavior is of the CFR Lorentz 
superdiffusive type for complete redistribution in the ARF but of the diffusive character if 
the excited atom is intrinsically coherent (even if the small scale resembles the pure Doppler 
superdiffusive character). So, one can use a generalized diffusive transport equation 5] for line 
radiation in this case, even if the effective mean transport coefficients will hide the details of 
the small scale partial redistribution effects. The present work is in this respect complemen- 
tary of the asymptotic analysis of the integral equation of line radiation transfer carried out 



by Frisch 
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33j. The work presented here is much simpler and straightforward than Frisch's 



approach. Furthermore, we clearly identify the CFR in the line core JW-type approximation 
for Rjj as incompatible with the large scale diffusive behavior. 

Figures [7] to [9] give the jump size PDFs computed numerically from Eq. (jlip . These give 
essentially the same information as the ODF PDFs, although in terms of the overall opacity 
length which is directly related to physical distances (for an additional note on the relation of 
opacity scales and physical distances, see ref. The superdiffusive behavior of R[ and Rm 

is evident, while the breakdown of super diffusion is manifested by the abrupt cut-off for Ru at 
higher jump lengths. Smaller details can also be pointed out. For Rj, there is an initial higher 
than superdiffusive Lorentz jump PDF which changes into the large scale Doppler asymptotic 
at high distances (transition not shown but confirmed for the last two cases of x' = 4 and 5 
in Fig. [7]). For the Rm redistribution, all cases converge into the same Lorentzian CFR. And 
finally, for the Rn case there are three characteristic jump length ranges; a first in which 
the jump PDF increases faster for decreasing opacities than CFR Lorentz, an intermediate 
with a CFR Doppler character (and these two distance ranges roughly correspond to the pure 
Doppler Ri case) and a third one corresponding to the diffusive onset due to a step cut-off 
for higher distances. 

The jump size PDFs define the overall opacity scale in which partial frequency redistri- 
bution effects are important as well as the transition from superdiffusion-like into standard 
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FIG. 7: Conditional one-sided Jump Probability Distribution Function (PDF) for Partial Redistribu- 
tion function i?j for several previous absorption frequencies x' in the sequence x 1 = 0, 1, 2, 3, 4 and 
5 (arrow direction). The Complete Redistribution limiting Doppler and Lorentz cases are also shown. 

diffusion for the Ru case. Prom Fig. [8] one expects that the onset of CFR Lorentz asymptotic 
for high opacities should roughly scale with the width a parameter; PFR effects should man- 
ifest themselves at distances up to 10 5 , 10 4 and 10 3 for a = 0.001, 0.01 and 0.1, respectively, 
when measured in a dimensionless overall opacity scale. On the other hand, in the case of Ru 
redistribution, Fig. [9] shows that the standard diffusion should become noticeable only at very 
high opacities; roug hly 10 6 , 10 5 and 10 4 for a = 0.001, 0.01 and 0.1, respectively. At smaller 
scales, and notably for finite systems with the biggest dimension not exceeding these opacity 
values, the faster than standard diffusion behavior is dominant. 

VIII. DISCUSSION AND CONCLUSIONS 

The line Opacity Distribution Function (ODF) approach was used to discuss the su- 
perdiffusive character of resonance line radiation transfer under Partial Frequency Redistri- 
butions (PFR) effects, as well as the conditions for its breakdown. These effects are the result 
of the ensemble average from the atom's rest frame (ARF) redistribution function into the 
lab frame (used to describe radiation transport). We discuss PFR effects for three frequency 
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FIG. 8: Conditional one-sided Jump Probability Distribution Function (PDF) for Partial Redistribu- 
tion function Rni for several previous absorption frequencies (x' = 0, 1, 2, 3, 4 and 5) and Voigt width 
parameters a. The Complete Redistribution limiting Doppler and Lorentz cases are also shown. 
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FIG. 9: Conditional one-sided Jump Probability Distribution Function (PDF) for Partial Redistribu- 
tion function Rn for several previous absorption frequencies x' and Voigt width parameters a. The 
arrow direction shown results in the sequence x' = 0, 1, 2, 3, 4 and 5. The Complete Redistribution 
limiting Doppler and Lorentz cases are also shown. 
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redistribution functions between absorption and reemission: Rj (pure Doppler broadening), 
Rji (Doppler plus natural) and Rni (combined Doppler, natural and collisional broadening). 
The Doppler ARF of Rj will give rise to the Complete Frequency Redistribution (CFR) su- 
perdiffusive Doppler asymptotic for sufficiently small line opacities, while the Lorentz ARF 
intrinsics of Rni give also small line opacity large scale Lorentzian CFR super diffusive asymp- 
totics. However, the ARF frequency coherence of Rn gives rise to a different behavior, be- 
cause the abrupt cut-off in the jump size distribution for higher distances means that there is 
a breakdown of the superdiffusive behavior. This is not recognized in both the astrophysical 
spectral line formation as well as in the lab scale atomic vapors studies. This breakdown 
means that the popular JW-type approximations of PFR under Rjj conditions are a severe 
approximation, since they assume that core radiation CFR allows a single scattering to pho- 
tons well into the absorption wings. Under experimental conditions rendering Rjj realistic, a 
generalized diffusion equation should be valid even for core line radiation. This has a classical 
analog found in physical kinetics, where the macroscopic hydrodynamic behavior is derived 
from a microscopic description through a Chapman-Enskog expansion l[. 

The Rjj and Rjjj redistribution functions assume either a complete frequency coherence or 
a complete redistribution in the ARF, respectively. These situations should correspond to two 
limiting cases for which the collisions during the excited state lifetime are almost nonexistent 
or very frequent. For resonance lines, the general form of the redistribution function should 



be 
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33, 



34] 



R (z'; x) = (1 - P c ) R n (x'; x) + P c Rjn (x'; x) , (40) 
where the branching ratio P c = r _pp — is defined as the probability that an elastic collision 

A C+-L rad 

will destroy the correlation between the absorbed and the reemitted frequencies in the ARF 
and should go to one in the limit of high densities and to zero in the low density limit. This 
equation reflects the competition between complete redistribution in the ARF on the one 
hand and almost complete coherence on the other hand, depending on the elastic collisional 
rate. Presumably, below some critical value for P c , Rjj-type redistribution will dictate a large 
scale diffusive behavior while after that critical value, Rjjj will take control of the large scale 
behavior giving rise to superdiffusion. The transition condensed in Eq. (I40p is important since 
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the classical landmark work of Post in the 1980s for trapping under PFR with 185 nm Hg 



radiation 



3l| corresponds to a range of P c values roughly from 1% to 30%. 



Figs. S] to [9] display the conditional ODFs and jump size distributions for a given previ- 
ous absorption optical frequency. The radiation migration in physical space is not a Markov 
process whenever PFR is meaningful since the whole hierarchy of the process cannot be 
constructed from the spatial distribution of the excitation in a given initial time plus the 
single (spatial) transition probability (contrary to the CFR case where the spatial transition 
probability is independent of past history). However, the spatial evolution of excitation can 
be embedded in a Markov process defined as the excitation migration in both spatial and 
frequency space [4J. This extended Markov process will consider explicitly the joint redistri- 
bution between absorption and reemission frequencies, information that otherwise would be 
contained implicitly in the past values of the absorption frequencies. In principle, it would 
even be possible to define two different Markov processes for the all the Rs redistribution 
mechanisms, depending on the coarsening level of the description [4]. On a fine level descrip- 
tion, one should have a complete stochastic formulation of the radiation trapping problem 
encompassing both the small scale (controlled by PFR) as well as the large scale (asymptotic) 
behaviors. But one could also devise a coarser level description retaining only the large scale 
asymptotic and blurring the small scale PFR details into a transition probability valid only 
for large distances. For the Ri and Rni cases, the large scale transition probability would 
be given by the corresponding Doppler or Lorentz CFR super diffusive asymptotic while for 
the Rn a simpler diffusion-type behavior should apply. 

Finally we discuss the applicability of the Opacity Distribution Function (ODF) formalism 
to astrophysical applications. In astrophysics, we have to deal with a very large number 
of spectral lines that overlap. Typically, consideration of all the lines (plus continua) in 
reaction-hydrodynamic modelling is prohibitive, but this can be circumvented by the use of 
an ODF. Detailed studies of this approach have shown that ODFs reproduce satisfactorily 
both emergent fluxes and the physical stellar atmospheric structure [26j]. In astrophysical 
radiative transfer problems, each individual opacity will be the result of grouping several 
lines. Recently, Wehrse and coworkers have used a Poisson distributed point process for 



the statistical description of the number of lines 



However, they did not attempt a 



comprehensive investigation of the dependencies of the ODFs on the detailed line profiles. 
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The present work is a contribution towards the improvement of this approach by including 
the effect of the absorption and reemission profiles in a similar way used by us to study the 
anomalous diffusion resulting from the redistribution along a single line profile. 
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